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Abstract 

A method is formulated for the description of lambda hypernuclei in the frame- 
work of the unitary-model-operator approach (UMOA). The method is applied to 
]\0. A lambda-nucleon effective interaction is derived, taking the coupling of the 
sigma-nucleon channel into account. The lambda single-particle energies are calcu- 
lated for the Osx/2) 0p 3 / 2 and 0p\/2 states employing the Nijmegen soft-core (NSC), 
Jiilich model-A (JA) and model-B (JB) hyperon-nucleon potentials. 
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1 Introduction 



There have been a number of attempts to study hypernuclei containing not 
only nucleons but also hyperons with strangeness [1-6] . It has been of high in- 
terest to investigate how the presence of strangeness broadens and/or modifies 
the knowledge achieved in the field of ordinary nuclear physics. 

It is a fundamental problem to understand properties of a quantum many- 
body system in terms of basic interactions between constituent particles. In 
hypernuclear physics a unified treatment of all the baryon-baryon interactions, 
especially nucleon-nucleon (NN) and hyperon-nucleon (YN) interactions, will 
be necessary if one wish to understand general properties of hypernuclei. How- 
ever, our knowledge on the YN interaction is quite inadequate in contrast to 
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that of the NN interaction due to experimental difficulties. In such a situa- 
tion many-body theoretical studies of hypernuclei could provide alternative 
information on the YN interaction. 

It has been discussed that physics of hypernuclei is different from conven- 
tional nuclear physics in some aspects. First, the channel couplings such as 
the EN- AN coupling play a significant role in the structure of hypernuclei in 
contrast to the very limited role of the AN-NN coupling appeared in nuclear 
many-body systems [7]. Second, the anti-symmetric spin-orbit force, which is 
forbidden in ordinary nuclei, emerges in hypernuclei in addition to the sym- 
metric spin-orbit force. 

The spin-orbit splitting of A single-particle levels in A hypernuclei was consid- 
ered experimentally to be much smaller than their nucleonic counterpart [8]. 
However, some recent studies as in Ref. [9] suggest that the splitting would not 
be so small as discussed earlier. It will be highly desirable for the spin-orbit 
splitting to be established in experiment. The magnitude of the spin-orbit 
splitting will put further constraints on the YN interaction. 

There have been several possibilities of solving the structure of hypernuclei. 
One possibility is to treat few-body or light hypernuclei by solving the Fad- 
deev equation [10] or by applying the cluster model [11,12]. Another possibility 
is microscopic studies of hypernuclei with a larger mass number. While a few 
lower partial waves play a decisive role in few-body or light hypernuclei, higher 
partial waves could give some contributions in heavier hypernuclei. These mi- 
croscopic studies of hypernuclei could make us understand interesting physics 
of their structure itself, and also help us in determining the YN interaction. 

Microscopic structure calculation of heavier hypernuclei is usually performed 
by introducing a YN effective interaction. For this purpose the G-matrix the- 
ory has been applied in two ways: In some works an effective YN interaction 
has been derived by simulating a nuclear matter G-matrix. In this approach 
the non-locality and energy dependence of the G-matrix are neglected, and 
the effective YN interaction is represented in the two- or three-range Gaus- 
sian form with the adjustable Fermi momentum kp. This type of the effective 
interaction has been referred to as the YNG interaction [13]. The YNG interac- 
tion has been successfully used in many of hypernuclear structure calculations 
[14]. It seems to us that the YNG interaction has made an important con- 
tribution to the considerable progress in overall description of hypernuclei. 
The other approach, in the framework of the G-matrix theory, is to calculate 
the G-matrix in finite hypernuclei by treating rigorously the Pauli exclusion 
operator and making a self-consistent calculation of single-particle potential 
of hyperons [15-19]. In contrast to studies of using the YNG interaction this 
approach possibly gives a better description of the state dependence of the 
effective YN interaction. It has been established that the G-matrix can be 
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used as a basic ingredient in constructing a microscopic many-body theory 
starting with bare interactions. 

In spite of a marked success in the applications of the G-matrix theory we 
may say from a formal point of view that the G-matrix is energy-dependent 
and non-hermitian, and does not have the property of decoupling between 
a model space of low- lying two-particle states and its complement. The G- 
matrix is not considered to be an effective interaction in a formal sense of 
the effective interaction theory. In order to derive an effective interaction we 
should add some higher-order corrections such as folded diagrams [20,21]. It 
would be desirable if we could have a theory describing many-body systems 
in terms of an energy-independent and hermitian effective interaction with 
the property of decoupling. Two of the present authors, K. S. and R. O., 
proposed a many-body theory, the unitary-model-operator approach (UMOA) 
that was formulated on the basis of such an effective interaction. The UMOA 
was applied to finite nuclei, 16 and 40 Ca [22,23], and some improvements have 
been attained. The method may be viewed as an alternative way of deriving 
an effective interaction in finite many-body systems. 

A unified study of hypernuclei and ordinary nuclei would be helpful in un- 
derstanding many-body systems of baryons in a microscopic way. Therefore, 
we extended the formulation of the UMOA, in the previous work [24], to the 
description of hypernuclei and made a calculation of the properties in l j0 with 
a realistic YN interaction. The purposes of this work are to present a general 
formulation of the extended UMOA and apply it to x JO by employing some 
realistic YN interactions. 

This paper is organized as follows: In Section 2 we present a general formula- 
tion of the UMOA for deriving the AN effective interaction. In Section 3 the 
approximation procedure is given for making actual calculations. In Section 4 
we apply the UMOA to x j0 by employing the Nijmegen soft-core (NSC) [25], 
Jiilich model-A (JA) and model-B (JB) [26] YN potentials. In Section 5 we 
make some concluding remarks. 



2 Formulation of the UMOA for the AN effective interaction 

We present a formulation of the UMOA for applying to a calculation of the 
properties in A hypernuclei, especially, closed-shell nucleus plus one A systems. 
A main purpose is to give a method of calculating an effective interaction 
between A and N . In the derivation of the AN effective interaction we should 
note that, because of the small mass difference of about 77MeV between A 
and S, the coupling of the AN and UN channels plays an important role. 
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We first consider a two-body hamiltonian of a UN- AN coupled system given 
by 



H YN = h YN + Am + v YN , 
where 



(2.1) 



hvM — 



o 



o 



(2.2) 



Am 



m E — m A 



and 



(2.3) 



Van-an Van-en 
Ven-an Ven-en 



(2.4) 



The terms tk and Uk for k = A, E and N are the kinetic and single-particle 
potential energies of A, S and N, respectively, and v Y n is the bare two-body 
YN interaction including the AN and EN channels. The terms m A and m E 
are the rest masses of A and E, respectively. The hamiltonian H YN describes a 
subsystem of interacting two particles AN or EN in the one-body potentials 
u A , u E and u N . 

In order to derive a AN effective interaction we introduce a unitarily trans- 
formed YN interaction given by 



v YN = e ~ SYN (h YN + Am + v YN )e SyN - (h YN + Am), (2.5) 

where S YN is the correlation operator defined in a space of two-particle states 
of AN and EN. The S YN is an anti-hermitian operator satisfying 



S YN — —S YN . (2-6) 

In the formulation of the UMOA the correlation operator S YN is determined 
by the equation of decoupling for v YN between a certain model space consisting 
of low- momentum two-particle states and its complement. This determination 
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of the correlation operator S YN is one of the characteristics of the present 
approach. 

For the determination of a AN effective interaction, we introduce a model 
space consisting of low-momentum AN states. Let P AN be the projection op- 
erator onto the model space. The complement of the P AN space includes a 
space of high-momentum AN states referred to as the Q AN space and a space 
of all the EN states referred to as the Q EN space. The condition of decoupling 
for the transformed interaction v YN is given by 



{Qan ~\~ Q en) v ynPan — PanV yn(Q AN ~\~ Q en) — 0. (2-7) 

Hereafter, we refer to the above equation as the decoupling equation. We 
assume that the one-body part h YN + Am in Eq. (2.5) is decoupled as 



(Qan "I" Q en ){h Y N "I" Am)P AN — Pan 

(h YN + Am) (Q AN + Q SN ) = 0. (2.8) 
On the above assumption the decoupling equation (2.7) becomes 

{Qan + QE N )e' SYN {h YN + Am + v YN )e s ™P AN = 0. (2.9) 

In general, the decoupling equation (2.9) does not determine S YN uniquely. 
The usual restrictive conditions are 

PanS yn P an = {Qan + Qen)S yn {Q an + Q en) = 0. (2.10) 

These conditions are called the minimal effect requirements [28], which mean 
that the transformation exp{S YN ) does not induce unnecessary transformation 
within each of the P AN and Q AN + Q EN spaces. 

In the general theory of constructing a hermitian effective interaction, the 
solution for S YN in Eq. (2.9) has been known [29-31] and is given by 

S YN = arctanh(w yiV — oo\ N ) 

{uj YN {u)i N uj YN ) n - h.c.} (2.11) 



OG 



v v 1)" 
^ 2n + l 



n=0 

with 



d 

U Y N = ^{QaN + Qen)\^ YNk)(^ANk\PAN, ( 2 -12) 
k=l 
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where d is the dimension of the P AN space. The \^ Y Nk) is the two-body eigen- 
state consisting of two components of AN and EN states, and satisfies 



YNk YNk I 



(2.13) 



where 



I* 



YNk) 



\^ANk) 
\^ENk) 



(2.14) 



The {i> AN k\ in Eq. (2.12) is the biorthogonal state of \ij) AN k) satisfying 



{4>ANk\4>ANk') — $kk', 
where \ip AN k) is the P^-space component of \^ Y Nk) defined by 



(2.15) 



ANk) 



Pan I ^ YNk) 



(2.16) 



It should be noted that the solution S YN in Eq- (2.11) is determined depen- 
dency on the choice of a set of d eigenstates {\^ YN k), k — 1, 2, d} as given 
in Eq. (2.12). We choose a set of {\^ Y Nk)} so that they have the largest P AN - 
space overlaps among all the eigenstates in Eq. (2.13). This is a usual choice 
in the derivation of the effective interaction. 

With the solution S YN in Eq. (2.11), the AN effective interaction in the P AN 
space is given by 



"Van — Pan^ynPa 



(2.17) 



It has been known that this AN effective interaction can be written explicitly 
[22,32] as 



\<pANa\v AN \<pANf3) 



an a I Pan \ yanP 

) + (l+^) 5 (0 

AN a \R\n\<^ AN0/ 



(1 (1+0?) 



(2.18) 



where \(p AN a) {\4>anp)) and ji a (fj,p) are defined through the eigenvalue equation 



^YN^YN \4>ANCt) — ^ a \4 > ANa)i 

and R AN is given by 



(2.19) 
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Ran — Pan(v Y n + V YN UJ YN )P AN . (2.20) 

The R AN is a P/iff-sp&ce operator that agrees with the effective interaction 
of non-hermitian type. Equation (2.18) is a formula for converting the non- 
hermitian effective interaction to the hermitian one. 

The term u A in Eq. (2.2) is introduced as a self-consistent potential of A 
calculated with the AN effective interaction v AN as 

(a A \u A \a' A ) = (^n\v an \o^ a ^ n ), (2.21) 

£/v:occupied 

where \a A ) (\a' A )) and |^) are the single-particle states of A and N, respec- 
tively. 

From Eqs. (2.18)-(2.20) we see that the AN effective interaction 

U AN Cctn DG 

calculated once the operator ui YN is given. The operator ui YN is determined 
completely, as given in Eq. (2.12), by solving the eigenvalue equation (2.13) 
and selecting d eigenstates {\^ YN k)} with the largest P^-space overlaps. It is 
noted here that Eqs. (2.11)-(2.21) give a set of equations that determines the 
AN effective interaction v AN self-consistently with the single-particle potential 
u A . 



3 Approximation procedure 

3.1 Two-step method for the calculation ofv YN 

Here we discuss a calculation procedure for obtaining the AN effective in- 
teraction acting in a small model space of low-lying AN states. We use the 
harmonic oscillator (h.o.) wave functions as basis states. It is desirable that 
the single-particle potentials of A, E and N in intermediate states are taken 
into account in a self-consistent way. However, as our main concern is to de- 
rive a AN effective interaction, the self-consistent treatment of u A would have 
predominant importance in the present calculation. For this reason we make 
a self-consistent calculation only for u A , and as for u N we employ a fixed po- 
tential calculated previously for closed shell nuclei in the UMOA [22]. The 
effect of the single-particle potential u E on the AN effective interaction has 
not necessarily been clarified, although some calculations with u E have been 
made for the G-matrices in hyperonic nuclear matter [33,34]. In the actual 
calculations made so far the single-particle potential u E has been considered 
to be rather shallow [14]. Although we have not known clearly how much the 
effect of u E on the AN effective interaction, we here assume u E = 0. 
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In solving the self-consistent equation for u A with v AN we should be careful in 
the following two respects: (a) In the calculation of energies of intermediate 
states in propagators we should use u A that is given for A states including 
high-momentum states. In order to do this we need a AN effective interaction 
acting in a large model space, because in general the Q AN space contains high- 
momentum AN states and the Q-space effective interaction Q AN v YN Q AN is not 
always well behaved. We thus make a calculation by introducing two model 
spaces, namely, large and small model spaces. The AN effective interaction in 
each of the model spaces is calculated by following a two-step procedure. In 
the first-step procedure, we define a large model space and calculate the AN 
effective interaction and single-particle potential of A in high-lying states. The 
AN effective interaction acting in this large model space will not be suitable for 
the effective interaction in a usual sense. By introducing a small model space 
and using the AN effective interaction determined in the first-step calculation, 
we proceed to the second-step calculation where we calculate a AN effective 
interaction acting in this small model space of low-lying AN states, (b) It is 
important in actual calculations how to determine the P and Q spaces because 
an effective interaction is determined dependently on the choice of the P and 
Q spaces. If a state in the P space mixes strongly with Q-space states in the 
eigenstate in Eq. (2.13), the matrix element of u YN or S YN defined in Sec. 2 
becomes large and as a result higher-order many-body correlations [35] give 
rise to non-negligible contributions. In general, the strong mixing takes place 
when some of the P- and Q-space states are quasi-degenerate in energy. For 
this reason, we should choose the P and Q spaces which are well separated in 
energy. 

With due regard to points (a) and (b) we perform a calculation of the effective 
interaction through the two-step procedure. In the following subsections we 
present the two-step method in detail. 



3.2 First-step decoupling 

Let us write a two-particle state of A and N as the product of the h.o. single- 
particle wave functions in the usual notation as 



Wy^n) = \n Y .l Y .j Y .m Y .,n N l N j N m N ), ( Yi = A, E). (3.1) 

We define a model space denoted by P^ and its complement Q A 1] N with a 
boundary number p 1 as 
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Px 



Pf 



o <1) 



N 



(b) 



Fig. 1. The model space and its complement in the first-step calculation. 
fp« if 2n A + l A + 2n N + l N < Pl , 

\olaPn) e < 

I otherwise. 
We also define a space of all the EN states denoted by Q { ^ N as 



(3.2) 



\ot s /3 N ) e Q 



(i) 

EN' 



(3.3) 



Note that the model space does not contain the UN states. 



If we make a graph of two-particle states, the and Q^ N spaces are 

represented in Fig. 1. The p x means a number such that (p 1 + 3)%{2/2 gives the 
maximum kinetic energy of the AN two-particle states in the Pj^ space. The 
JiQ denotes the unit of the h.o. energy which we use commonly for A, U and 
N. As we require to introduce a sufficiently large model space in the first-step 
calculation, we should choose p x as large as possible. Another requirement is 
that the states in the P^ space should be separated in energy from the states 
in the Q A 1] N space. For this reason, we adopt a model space Pj^ of a triangular 
shape. This choice of the Pj^ space guarantees the separation in energy. 

In Fig. 1 the spaces Q^ N for Y { = A and U, specified by the numbers p 1 ,p F 
and p x , stand for the excluded spaces due to the Pauli principle for nucleons. 
The number p F is defined as p F = 2n N + l N with the quantum numbers {n N , l N } 
of the highest occupied orbit in the closed-shell core nucleus. The p x means a 
number such that (p x + 3)hf2/2 gives the maximum kinetic energy of the YN 
two-particle states in the and spaces. We then write the Qffi and 
spaces as 



\a A (3 N ) e Q 



(X) 
AN 



if p 1 < 2n A + l A + 2n N + l N < p x 
and < 2n N + l N < p F 



(3.4) 
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and 



\u s f3 N ) G if < 2n s + l E + 2n N + l N < p x 

and < 2n N + l N < p F . (3.5) 

We also define the sum of the excluded spaces as 



Q% = Q% + Q%- (3-6) 

Furthermore, by using the projection operators, Q^l, Q ( ani Qen an d Q ( eni we 
introduce projection operators onto the spaces of allowed YN states as 



Qan — Qan Qani (3-7) 
Qsn = Qsn ~ Qeni (3-8) 

and we also define 

Qyn = Qan + Qen- (3-9) 

Considering the Pauli-exclusion principle for nucleons the terms Q^v Y nQ { yn 
and P^v Y nQ ( yn + n - c - should be removed in treating the problem of two-body 
correlations. However, if we remove all the matrix elements of these terms, this 
leads to overmuch counting of the Pauli principle. The matrix elements which 
are diagonal in nucleon states should be restored even though nucleons are 
in occupied states. Therefore, instead of v YN , we should use a YN interaction 
v YN defined as 

v yn ^(p^+q ( z>y N (p^+q? n )+y: e q?»w%q>). 

Vn Y i= A,E 
Y] =A,E 

(3.10) 

where v N denotes the set of the h.o. quantum numbers {n N , l N } of an occupied 
state of a nucleon, and the operator Q^ N (v N ) is defined as 

Q ( ^K) = EK^)K^I' (Y t = A,E). (3.11) 

Oty- 
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The two-body YN equation to be solved in the first-step calculation is now 
written as 



+ Am + v YN )\ YN;h) = E kl \YN;h), (3.12) 



where we define h^l as 



AS _ + < + + | (3.13) 

t E + t N + u N 

The u A 1] is the single-particle potential of A to be calculated with the AN effec- 
tive interaction self-consistently. The number k\ stands for a set of quantum 
numbers to specify a two-body YN eigenstate. The unitarily transformed YN 
interaction in the first-step calculation is given by 



vW = e^(/i« + Am + v YN )e s ™ - (h% + Am). (3.14) 
The decoupling equation for determining the correlation operator S ( yl becomes 

Qyn^yn-Pan = Pan^ynQyn = 0- (3.15) 
The AN effective interaction v^l in the P$ space is then given by 

< = (3-16) 

In the first-step procedure we solve the decoupling equation (3.15) by intro- 
ducing some assumptions. As we take a sufficiently large number as p 1; we 
may assume that the one-body potential of A for a state in the Q% space is 
negligible in comparison with the sum of the energies, t A + t N + u N . Therefore, 
takes nonzero values only for states in the Pfy space, and we have 

= PXPZ (3.17) 

The definitions of Pj£> and Q^ N in Eq. (3.2) allow us to divide a space of AN 
states written in relative and center-of-mass (cm.) coordinates into the same 
subspaces P^ and Q^ N as 



\AN;nlSJ r ,N c L c ) e < 



P« ii2n + l + 2N e + L e <p u 

(3.18) 

Qy N otherwise, 
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where {n, 1} and {N c , L c } are the h.o. quantum numbers of relative and cm. 
motions, respectively. The S and J r denote, respectively, spin and angular 
momentum defined through J r = I + S. The isospin T is omitted because T 
takes a constant 1/2 in the AN system. A EN state in the Q^ N space in Fig. 1 
is also written by the product of the h.o. relative and cm. states. 

The YN interaction v yn in Eq. (3.10) is not diagonal in the cm. quantum 
numbers {N c , L c } because the projection operators Pj^, Q YN and Q ( Y . N are n °t 
diagonal in {N C ,L C }. Thus, the AN effective interaction vfy is not diagonal 
in {N C ,L C }. It has been shown, however, in the calculation of ordinary nuclei 
that {N C ,L C } non-diagonal elements of the G-matrix are generally small [36]. 
Therefore, we make an assumption that the {N c , L c } non-diagonal matrix 
elements of vfy can be neglected. 

The two-body eigenvalue equation (3.12) is solved for the relative motion with 
a fixed set of {N C ,L C } written by 



(h^ + Am + v YN )\ YN; k u ll'SS'J r , N C L C ) = E kl | YN; k u ll'SS'J r , N C L C ), 

(3.19) 

where we rewrite h^l as 



h (1) 



t r AN ~\~ t c AN ~\~ U A ^ + U 



YN 



t r SN 4" t c EN ~\~ U N J 



(3.20) 



TllG t r YN cLIld t c YN 

are the kinetic energies of the YN relative and cm. motions, 
respectively. The state | YN; k±, ll'SS'J r , N C L C ) expresses a superposition of the 
AN and EN h.o. states and gives a general form of the eigenstate when we 
consider the EN-AN coupling and the YN interaction with the tensor and 
antisymmetric spin-orbit couplings. The ki is an additional quantum number 
specifying a two-body YN eigenstate. 

In general, the one-body part + Am does not satisfies the assumption of 
decoupling in Eq. (2.8). However, it has been shown in the calculation in nuclei 
that if we take a larger number as p 1 , the effect of violating the assumption of 
decoupling becomes considerably smaller [22,23]. We have confirmed that this 
non-decoupling effect becomes negligible in calculating nuclear effective inter- 
actions. We here assume this to be the case also in the determination of the 
AN effective interaction. Furthermore, we introduce another approximation 
that the sum of t cY%N and u Y . + u N be diagonal in the cm. quantum numbers 
N c and L c . This assumption has also been made sure to be acceptable in the 
previous works [22-24]. 

In order to solve the eigenvalue equation (3.19) we need to calculate the ma- 
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trix elements of uty + u N in the relative and cm. coordinates representation. 
However, the calculation of these matrix elements is difficult because the term 
iiy' + u N can not be separated into a sum of operators written in relative 
and cm. coordinates. Therefore, as have made in the previous work [22], we 
employ the angle-average approximation and have 



+{n'l'N c L c : X\n Y .l Y .n' N l N : X) YiN {n N l N \u N \n'J N ) )5 W 5 SS ', (Y { = A, E), 



where, as was assumed previously, uty = for Y, = S. The {nlN c L c : 
\\n Y .l Y .n N l N : \) YiN is the h.o. transformation bracket [37] for the YN two- 
body state which depends on the rest masses for Yj = A and S. Here the 
angle- average approximation means that, when a function f(jm) is given with 
an angular momentum j and its ^-component to, we approximate f(jm) as 
f(j) — Em /(j m )/(2j + 1)- In the present case the matrix element on the 
left-hand side in Eq. (3.21) is dependent on the ^-components m Jr and m Lc of 
J r and L c , respectively. The average over m Jr and m Lc leads to the expression 
given in Eq. (3.21). This approximation is essentially the same as has been 
made to derive an angle-averaged Pauli operator by Wong and Sauer [38,39]. 

The use of the basis states {| FjiV; nlSJ r , N C L C )} makes it difficult to represent 
the operator Q ( YN , but it has been confirmed that the angle-average technique 
works well. We thus make an approximation to the operator Q ( y N as 




(3.21) 




£ 9 YiN (nl, N C L C )\ YN; nlSJ r , N c L c )(Y t N; nlSJ r , N C L, 



nlN c L c SJ r 
Yi=A,S 



(3.22) 



where the summation should be made on the conditions 



p 1 <2n + l + 2N C + L c < p x for Y t = A, 
< In + I + 2N C + L c < p x for Y, L = S. 



(3.23) 



Here the coefficient 9 Y . N (nl, N C L C ) is defined as 



9 Y . N (nl,N c L c ) 
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= 1 - E 7WTT^T~rT\ ( nlNcLc : X \n Y l Yi n N l N : A) 2 



(Y l = A,E), (3.24) 
where in the summation the following conditions should be satisfied 

< 2n N + l N < p F , 

p 1 < 2n Y . + l Y . + 2n N + l N < p x for Y, = /l, 

< 2n Y% + Z y . + 2n w + l N < p x for Y" 8 = £. (3.25) 

Furthermore, we need an expression of Q Y l(v N ) given in Eq. (3.11) in the 
angle-average approximation 



- E e' Y . N {nl,N c L c ,nJ N )\Y i N-nlSJ r ,N c L c ){Y i N-nlSJ r ,N c L c l 

nlN c L c SJ r 

(Yi = A,Z), (3.26) 

where the summation should be made on the same conditions as in Eq. (3.23). 
The coefficient 9 YN (nl, N c L c ,n N l N ) is given by 



0' Y , N (nl,N c L c ,n N l N ) 

= E / , ,!to/ ,.A nlN c L c : X\n Y .l Y .n N l N : A) 2 ( = /l, 17), 

(3.27) 

where n Y% and Z y . should satisfy the conditions given in Eq. (3.25). 

Once the set of the eigenstates | YN; ki, U'SS'J r , N C L C ) in Eq. (3.19) is solved, 
the AN effective interaction acting in the Pj^ space is given according to 
Eqs. (2.18)-(2.20) by substituting uj yn = w^, where oo { Y ] N is defined as 



" { ?» = E Qyn\ YN; ki, U'SS'Jr, N C L C )(YN; h, WSS'J r , N C L C \P™. 

ki,ll'SS'J r ,N c L c 

(3.28) 

With the AN effective interaction v A 1 ^ thus determined, we calculate a new 
single-particle potential {n A l A \u A 1) \n' A l A ) as an averaged matrix element of {(n A 
IaJaW^WaJa), 3a = h ± 1/2}. The (n A l A \u^\n'J A ) can be given explicitly in 
terms of the effective interaction v^l as 
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\n A t, A \u A \n A L A ) — — — - — - 

^=^±1(2^ + 1) 

_ (2A + l)(2J r + 1) 

" nnnsh^ (2/, + l)(2/ + l) 

{n N j'jv } :occupicd 

x (nlN c L c : \\n A l A n N l N : X) AN (n'lN c L c : X\n'J A n N l N : A)^ 
x (AN; nlSJ r , N c L c \d A 1} N \AN; n'lSJ r , N C L C ). (3.29) 

In the above treatment of the single-particle potential of A in the first-step 
calculation, we do not consider the one-body spin-orbit splitting. 

The AN effective interaction v A 1] N in Eq. (3.16) is dependent on the single- 
particle potential u A l \ and u A r> is determined with v^l as in Eq. (3.29). There- 
fore, the equation for determining v A 1] N and u A 1} should be solved simultaneously. 



3.3 Second-step decoupling 



The model space Pj^ in the first-step procedure is chosen as a low-momentum 
space, but the Pj^ space is still large when we wish to use v ( A l in the calculation 
of bound or low excited states of A hypernuclei. With v^l, we further proceed 
to the calculation of the AN effective interaction acting in a smaller model 
space and correspondingly the single-particle potential of A in low-lying states. 
We separate the Pj^ space into the Pj^, Pj*\ Pj$ and Q ( A } N spaces specified 
by the numbers p u p F , and p 2 as shown in Fig. 2. The number p 2 is introduced 
to specify the uppermost bound state of A. Here these spaces are defined as 



\a A p A 



e < 



P™ if 2n A + l A < p 2 and 2n N + l N < p 



F- 



Pjy if 2n A + l A > p 2 , 2n N + l N < p F and 2n A + l A + 2n N + l N < p x , 
P ( AN ] if 2n A + l A < p 2 , 2n N + l N > p F and 2n A + l A + 2n N + l N < p 1 , 
Q A 2) N if 2n A + l A > p 2 , 2n N + l N > p F and 2n A + l A + 2n N + l N < p 1 . 

(3.30) 



The Pjy space is the excluded space due to the Pauli principle for nucle- 
ons. We are principally interested in the AN effective interaction that gives a 
description of the bound states of A. It would be desirable that such a AN ef- 
fective interaction satisfies the condition of decoupling between two AN states 
consisting of the bound and unbound states of A. This requirement is equiv- 
alent to the situation that the AN effective interactions should be decoupled 
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A 



Pi 



P 2 



p (2) 
TAN 



n (2) 

Wan 



p< Y > 
Tan 



Pi 



N 



Fig. 2. The model space and its complement in the second-step calculation. 

between two AN states in the Pj® and Q (2) N spaces. Therefore, we solve the 
decoupling equation for v {2) N between the AN states in the Pf^ and Q (2) N spaces. 

The space Pj$ is not a space to be excluded due to the Pauli principle. How- 
ever, some states in the Pj$ and Pj^ spaces overlap each other in energy. 
Therefore, we treat separately the effect of the coupling between the P}$ and 
Pj^ spaces. Actually the interaction between states in the P$ and P {2 ^ spaces 
induces the core polarization in the core nucleus, as discussed in Refs. [22] and 
[24] . We take into account the effect of the core polarization in the usual per- 
turbative calculation. 

In the second-step calculation the unitarily transformed interaction is given 
by 



-,(2) 



;(2) 



+ v (1) 

AN ' AN 



;(2) 



- h (2) 

AN ' 



(3.31) 



where 



h (2) 

AN 



t A + u (2) +t N + U N . 



(3.32) 



The term u A 2) is the self-consistent potential of A to be determined in the 
second-step procedure. The decoupling equation for determining the correla- 
tion operator S (2 l in Eq. (3.31) becomes 



r)( 2 ) 7~>< 2 ) P< 2 ) — P(2)f,(2) /Q< 2 > _ f| 
Ka CAN u AN ± AN ± AN u AN Ka CAN u 



(3.33) 



We note here that the assumption of decoupling 



Q 



(2) 1,(2) D(2) _ r. 
AN^AN 1 AN u 



(3.34) 
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is satisfied exactly because the two two-body states in the P}$ and Q A 2) N spaces 
do not contain a common single-particle state of A or N. 

The two-body eigenvalue equation to be solved in the second-step calculation 
is written as 



(^ + Q%)(h% + i;<£)(P$ + Qf N )\AN; k 2 , J) = E k2 \AN; k 2 , J), (3.35) 

where J denotes the total angular momentum, and k 2 a set of quantum num- 
bers of a AN two-body eigenstate. It should be noted that this eigenvalue 
equation can be solved exactly in the P^ + Q (2) N space by the diagonalization 
of the hamiltonian written with the basis states {\n A l A j A ,n N l N j N ; JM)} if the 
Pan + Qan space is chosen suitably as a small space. Once a set of the eigen- 
states \AN;k 2 , J) in Eq. (3.35) is solved, the AN effective interaction v A 2 l in 
the P {2 ^ space is given from Eqs. (2.18)-(2.20) with oj yn = uj a 2) n defined by 



< = E Qa 2) n\AN; k 2 , J)(AN; k 2 , J\P%. (3.36) 

k2,.J 



The single-particle potential u A 2) is calculated with the AN effective interac- 
tions, v ( AN and v ( AN , as 

(a A \uf\a' A ) 

{a A (3 N \vf N \a' A (5 N ) if \a A (3 N ) e P { AN and \a' A f3 N ) G 



^ E 

Pn :occupied 



(a A /3 n\^ < an\ c< 'aPn) otherwise. 



(3.37) 



We solve the set of equations for determining v {2) N and u A 2) iteratively until the 
result converges. 

We remark here that we employ v ( A l as the effective interaction in the Pj^ 
space in the calculation of u (2> . Since the effective interaction in the Pj$ space 
is unchanged in the transformation made in the second-step procedure, we 
may write the AN effective interaction in the space as P^v^P^ = 
Pan^anPan- We, therefore, have used partly the AN effective interaction v ( AN 
in the calculation of u A 2) as in Eq. (3.37). 

3.4 Corrections in perturbative treatment 



First we note that the one-body hamiltonian h A 2) N in Eq. (3.32) has non- 
diagonal terms because we did not make the Hartree-Fock calculation. These 
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Fig. 3. Calculated diagrams for the A single-particle energies. The AN effective inter- 
action v AN is represented by the wavy-line vertex, and the x vertex with the dashed 
line represents the non-diagonal part of the one-body hamiltonian = + 
for k = A and N. The propagations of N and A are represented by the single and 
double external lines, respectively. 

terms could bring about sizable corrections to the single-particle energies of 
A. We adopt the same h.o. energy both for A and N, though it has been 
considered that the wave function of A is basically different from that of N. 
Therefore, we must evaluate the corrections arising from the non-diagonal 
terms of h A 2) N . By taking the diagonal part of h^l as the unperturbed hamilto- 
nian we calculate these corrections perturbatively up to second order as shown 
in diagram (a) of Fig. 3. 

Next, the interaction between the AN states in the P^ and P^Q spaces in- 
duces the core polarization in the core nucleus, and is important in the cal- 
culation of the corrections. We here assume that these corrections can be 
calculated approximately, using the interaction Pj^v^P^ + h.c. It will be 
desirable to take into account the effect of intermediate nucleon states with 
high momentum. In this perturbative calculation we include only the nucleon 
states belonging to the P$ space, that is, the nucleon h.o. state should satisfy 
2n N + l N < Pi — (2n A + l A ), where {n A ,l A } are the h.o. quantum numbers of 
a A state. The validity of this truncation of nucleon states will be checked in 
the numerical calculation in the next section by increasing the number p 1 . The 
diagrams corresponding these corrections are given in diagram (b) and (c) of 
Fig. 3. 



4 Application to 



We performed calculations employing the Julich model-A (JA), model-B (JB) 
[26] Q and Nijmegen soft-core (NSC) [25] [^potentials for the YN interaction. 

1 We used the computer code for the Julich YN potentials in the momentum-space 
representation made by B. Holzenkamp and A. Reuber. 

2 We used the computer code for the Nijmegen soft-core YN potential in the 
coordinate-space representation provided by T. A. Rijken and J. J. de Swart. 
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Fig. 4. The p Y dependence of the calculated A single-particle energies including the 
perturbation corrections given in Fig. 3. 

We used the h.o. basis states with a common UQ = 14MeV for A, E and N. 
The dependence of the effective interaction and single-particle energies on the 
value UQ has been already discussed in 16 O [35] and we confirmed that the 
dependence was quite small around %Q = 14MeV if perturbation corrections 
were included. As for the nucleon single-particle potential u N , we used the 
fixed data [22] calculated using the Paris potential [27], and we have assumed 
that u E = 0. The boundary numbers p 2 , p F and p x are taken as p 2 — 1, p F — 1 
and p x = 12 for l JO. 

We present, in Fig. 4, the p 1 dependence of the single-particle energies of 
A including the perturbation corrections shown in Fig. 3. The parameter p 1 
specifies the Pjfy space in which the AN effective interaction vfl is determined 
rigorously without any approximation as discussed in Sec. 3. Therefore, as p 1 
becomes larger, the calculated result does more accurate. If the result con- 
verges at a certain number p u we may say that the method of the two-step 
calculation of the AN effective interaction would be trustworthy. As shown in 
Fig. 4, the results are convergent and stable for p 1 > 7 for each YN potential. 
If we take p 1 larger than 6, our model space is considered to be sufficiently 
large. We note that the model space with p 1 — 6 corresponds to that used 
in the G-matrix calculations made by Halderson [16] and by Hao and Kuo 
[17,21]. 

Table 1 shows the dependence of typical matrix elements of the AN effective 
interaction at p x — 8 on the value N max which denotes the number of the 
h.o. basis functions. The number iV max determines the dimension of the space 
of relative states in the first-step calculation. When we solve the eigenvalue 
equation (3.19), taking into account the tensor and UN- AN couplings, we need 
2x2xiV max basis functions in each partial- wave channel. Table 1 shows that 
the results are convergent if we take 120 for the JA and JB potentials, and 
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200 for the NSC potential as N max . We have also made sure that the situation 
of convergence is similar in the other matrix elements of the AN effective 
interaction. As we employ the NSC potential given in the configuration-space 
representation, we need much more basis states than in the Jiilich potential 
given in the momentum-space representation in order to treat short-range 
correlations accurately. 

Table 1 

Dependence of typical diagonal matrix elements (ab\v^l\ab) j on A max in units of 
MeV, where a and b denote the A and N single-particle orbits, respectively, labeled 
as 1 = 0si/2) 2 = Opi/2 arid 3 = OP3/2, and J is the total angular momentum. 





a = 2, 


b = 2, 


J= 1 


a = 3, 


6 = 2, 


J = 2 


N max 


JA 


JB 


NSC 


JA 


JB 


NSC 


40 


-0.991 


-1.331 


-0.347 


-0.825 


-0.828 


-1.146 


60 


-1.100 


-1.432 


-0.402 


-0.909 


-0.892 


-1.186 


80 


-1.145 


-1.483 


-0.408 


-0.932 


-0.920 


-1.191 


100 


-1.152 


-1.492 


-0.412 


-0.937 


-0.926 


-1.195 


120 


-1.152 


-1.492 


-0.417 


-0.937 


-0.925 


-1.198 


140 






-0.422 






-1.201 


160 






-0.426 






-1.204 


180 






-0.428 






-1.205 


200 






-0.428 






-1.205 



Table 2 

Diagonal matrix elements (ab|i5^|a6)j in units of MeV. The notations are the same 
in TableL 



a 


b 


J 


JA 


JB 


NSC 


a 


b 


J 


JA 


JB 


NSC 


1 


1 





-1.31 


-0.06 


-4.04 


2 


3 


1 


-0.38 


-0.75 


-0.42 


1 


1 


1 


-2.49 


-3.47 


-1.38 


2 


3 


2 


-0.94 


-0.92 


-1.26 


1 


2 





-0.78 


-1.22 


-0.39 


3 


1 


1 


-0.94 


-0.62 


-1.45 


1 


2 


1 


-0.95 


-0.98 


-0.70 


3 


1 


2 


-1.38 


-1.80 


-1.06 


1 


3 


1 


-0.46 


-0.11 


-1.13 


3 


2 


1 


-0.50 


-0.90 


-0.34 


1 


3 


2 


-1.17 


-1.50 


-1.00 


3 


2 


2 


-0.94 


-0.93 


-1.21 


2 


1 





-1.05 


-1.56 


-0.54 


3 


3 





-0.49 


0.65 


-2.03 


2 


1 


1 


-0.84 


-0.84 


-0.95 


3 


3 


1 


-0.73 


-1.16 


-0.02 


2 


2 





-0.10 


0.57 


-0.71 


3 


3 


2 


-0.20 


0.12 


-0.96 


2 


2 


1 


-1.15 


-1.49 


-0.43 


3 


3 


3 


-1.34 


-1.78 


-0.84 
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We show, in Table 2, the diagonal matrix elements of the Os-Op shell AN effec- 
tive interaction calculated in the second-step procedure. Because of the differ- 
ent characteristics of the original YN potentials, the results depend strongly 
on the YN potentials used in the calculation. Considering the situations of 
convergence shown in Fig. 4 and Table 1, we have used the numbers p x = 8, 
and iV max = 120 and 200 for the Jiilich and Nijmegen potentials, respectively, 
in Table 2. We shall use hereafter the same set of the numbers in the following 
tables. 

In the present calculation we consider three kinds of the treatments of the 
Pauli-blo eking effect in the first-step procedure. 

(1) The Pauli blocking in the spaces Q<$ and is not considered, that 
is, the original YN interaction v YN is used in place of the interaction v YN in 
Eq. (3.10). 

(2) The Pauli blocking in both the spaces and is considered, that is, 
we take only the first term of v YN in Eq. (3.10). 

(3) The Pauli blocking in the spaces and is considered, but the 
diagonal interaction, the second term of v YN in Eq. (3.10), is restored. 

It is noted that, in the above three cases, the predominant Pauli-blocking 
effect which comes from the low-momentum AN states contained in the Pj^ 
space is taken into account rigorously in the second-step calculation. 

We here discuss the dependence of the different treatments of the Pauli- 
blocking effect in the first-step procedure for three YN interactions. In Table 3 
we display the single-particle energies of A including the perturbation correc- 
tions shown in Fig. 3. Rows (1), (2) and (3) correspond to different treatments 
(1), (2) and (3) of the Pauli-blocking effect in the first-step procedure. All the 
single-particle energies of A in the 0si/ 2 , 0p 3 / 2 and 0pi/ 2 states in row (1) are 
negative for three kinds of the YN interactions, and the results for the JB 
and NSC potentials have almost the same order of magnitude. The results 

Table 3 

Dependence of the single-particle energies of A on the different treatments of the 
Pauli-blocking effects for the Osi/2, 0^3/2 and Opi/2 states including the perturbation 
corrections. 







OS1/2 






0p 3 /2 






OPl/2 




Pauli 


JA 


JB 


NSC 


JA 


JB 


NSC 


JA 


JB 


NSC 


(1) 


-14.09 


-16.47 


-16.74 


-2.49 


-3.95 


-4.47 


-1.29 


-2.46 


-3.07 


(2) 


-10.96 


-13.42 


-9.23 


-0.16 


-1.76 


-0.04 


1.28 


-0.17 


1.68 


(3) 


-12.19 


-15.54 


-11.69 


-0.78 


-2.81 


-1.17 


0.64 


-1.24 


0.49 
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in row (2) show that all the single-particle energies are less attractive than 
those given in row (1), and especially, those for the Opi/2 state are positive 
(unbound) for the JA and NSC potentials. This means that the Pauli-blocking 
effect contributes repulsively to the single-particle energies and plays a very 
significant role although it depends on the YN potential used. In row (3) the 
single-particle energies become attractive in comparison with those given in 
row (2) because of the restoration of the "diagonal" interaction. The results 
in row (3) for the JA and NSC potentials have almost the same order of mag- 
nitude in contrast with those given in row (1). On the whole, we see from 
Table 3 that the Pauli-blocking effect is significantly large and yields different 
contributions to the single-particle energies dependently on the YN potentials. 
The difference between the results shown in rows (2) and (3) is fairly large. 
Therefore, we must treat the Pauli-blocking effect carefully and should adopt 
the treatment of case (3) as the most accurate one. In Fig. 4, Table 1 and 
Table 2, we have displayed the results of the treatment of the Pauli-blocking 
effect in case (3). 

Next, the partial- wave contributions to the first-order potential energy of A 
in the Osi/2 state are given for the different treatments of the Pauli-blocking 
effect. We see from Table 4 that the contributions of the 1 S , and 3 Si- 3 D 1 
channels are very different dependently on the YN interactions employed. In 
the 1 Sq channel the largest Pauli-blocking effect is brought about in the result 
for the JA potential. On the other hand, in the 3 Si- 3 Di channel we have an 
extremely large Pauli-blocking effect for the NSC potential. This result shows 
that the NSC potential has the strong EN- AN coupling in the 3 Si- 3 Di channel 
as discussed in Ref. 20). The absolute values contributing to the first-order 
potential energy of A in the channels with angular momenta I > 1 are much 

Table 4 

Partial wave contributions to the first-order potential energy of A in the 0s 1 / 2 state. 
All entries are in MeV. 









3 S 1 - 3 D 1 


3 Po 


l Pi- 3 Pi 


3 P 2 - 3 F 2 


Total 




(1) 


-4.87 


-17.94 


0.32 


0.91 


-0.10 


-21.67 


JA 


(2) 


-2.51 


-17.26 


0.40 


1.14 


-0.09 


-18.32 




(3) 


-3.14 


-17.82 


0.37 


1.07 


-0.10 


-19.61 




(1) 


-0.24 


-24.81 


0.31 


1.43 


0.31 


-23.01 


JB 


(2) 


0.12 


-22.61 


0.37 


1.59 


0.32 


-20.22 




(3) 


-0.17 


-24.43 


0.35 


1.54 


0.31 


-22.40 




(1) 


-9.70 


-14.16 


0.18 


1.80 


-1.86 


-23.73 


NSC 


(2) 


-9.58 


-7.29 


0.28 


2.05 


-1.84 


-16.38 




(3) 


-9.59 


-9.79 


0.24 


1.96 


-1.85 


-19.04 
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smaller than those in the 1 S and 3 Si- 3 Di channels. 

In Table 5 we display the single-particle energies of A with the contributions 
of the perturbation corrections given in Fig. 3, where the Pauli-blo eking effect 
is treated accordingly to case (3). From this table we see that the J A and NSC 

Table 5 

Contributions to the single-particle energies of A in units of MeV. Here KE and 
PE stand for the kinetic energy and the first-order potential energy, respectively. 
Rows (a), (b) and (c) are the contributions of diagrams (a), (b) and (c) in Fig. 2, 
respectively. The quantity Ag; s is the spin-orbit splitting of the Op states. 







JA 








0*1/2 


0^3/2 


Opi/2 


As/, 


KE 


10.50 


17.50 


17.50 




PE 


-19.62 


-14.50 


-12.97 


1.53 


(a) 


-1.73 


-2.83 


-3.03 


-0.20 


(b) 


-1.00 


-0.43 


-0.05 


0.39 


(c) 


-0.33 


-0.51 


-0.81 


-0.31 


Total 


-12.19 


-0.78 


0.64 


1.41 






JB 








Osi/2 




Opi/2 


Ae, s 


KE 


10.50 


17.50 


17.50 




PE 


-22.44 


-16.43 


-14.81 


1.62 


(a) 


-1.39 


-2.19 


-2.36 


-0.17 


(b) 


-1.31 


-0.51 


-0.11 


0.40 


(c) 


-0.90 


-1.17 


-1.45 


-0.28 


Total 


-15.54 


-2.81 


-1.24 


1.57 






NSC 








OSl/2 


0^3/2 


QPl/2 


Ae* s 


KE 


10.50 


17.50 


17.50 




PE 


-19.06 


-14.71 


-12.95 


1.76 


(a) 


-1.68 


-2.76 


-3.04 


-0.28 


(b) 


-0.93 


-0.52 


-0.12 


0.41 


(c) 


-0.52 


-0.68 


-0.90 


-0.22 


Total 


-11.69 


-1.17 


0.49 


1.66 
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potentials yield almost the same values of the first-order potential energy PE 
and the perturbation corrections in spite of having different partial-wave con- 
tributions as shown in Table 4. The single-particle energies obtained for the 
JB potential are most attractive among the YN interactions used. It should be 
noted here that the values of the KE plus PE in the Op states for three kinds of 
YN interactions are positive (unbound). Therefore, the perturbation correc- 
tions play a decisive role when we argue that the A in the Op states is bound 
or not. Correction (a) in Fig. 3 is induced by the presence of non-diagonal 
terms in the one-body hamiltonian h A — t A + u A . This effect should be treated 
in the Hartree-Fock (HF) calculation. However, in the present calculation, we 
treat this effect perturbatively. These non-diagonal terms work to change the 
single-particle wave function of A. Therefore, this procedure is indispensable 
in this approach since we have taken a fixed value of UQ commonly for A 
and N. This correction (a) contributes most attractively to the single-particle 
energies among all the corrections in Fig. 3. The magnitude of every correc- 
tion is considerably small in comparison with the first-order potential energies 
for all the YN interactions used. All the perturbation corrections contribute, 
however, attractively to the single-particle energies of A, and do not change 
the mutual spacings of the single-particle states significantly. This trend is 
very different from the situation observed in ordinary nuclei such as 16 O and 
40 Ca, in which the correction terms contribute to about 40% of the calculated 
spin-orbit splitting of the Op or Od states [22,23]. The role of many-body cor- 
relations, such as the core polarization, does not seem to be important in the 
spin-orbit splitting of A. 



5 Concluding remarks 

We have presented a formulation of the UMOA to apply it to a structure cal- 
culation of hypernuclei. We have introduced a unitary transformation exp(S'), 
with the correlation operator S, that describes two-body correlation of hy- 
peron and nucleon. An equation has been given for determining the correla- 
tion operator S as the equation of decoupling so that the unitarily transformed 
interaction does not have non-zero matrix elements between the model space 
P of low-lying AN states and its complement Q. The effective interaction has 
been given by the projected interaction onto the P space. The effective inter- 
action thus determined has properties of being hermitian, energy-independent 
and decoupled between the P and Q spaces. In comparison with the usual 
G-matrix approach, the use of such an effective interaction would have some 
advantages in describing many-body systems. 

The UMOA has been applied to the calculation of the structure of ]\0. The 
AN effective interaction has been calculated by employing three YN poten- 
tials, namely, the JA, JB and NSC potentials. With the AN effective inter- 
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action, we have calculated the single-particle energies of A in the 0si/ 2 , 0p 3 / 2 
and Opi/2 states. 

We have seen that how to treat the Pauli-blocking effect is very important, 
especially when we solve the problem of the EN-AN coupling. In the usual 
treatment one excludes all the interactions acting in a space of AN and EN 
states in which a nucleon is in occupied states. However, this treatment of the 
Pauli-blocking effect would lead to overmuch counting. We have introduced 
a new treatment in which, while we exclude the interactions acting among 
forbidden YN states as in the usual treatment, we restore the interactions that 
are diagonal in occupied nucleon states. We have observed that this treatment 
of the Pauli-blocking effect has given rise to a significant effect as compared 
with the usual treatment. 

The present results have shown that A in the Opi/2 state is unbound for the J A 
and NSC potentials and, on the other hand, bound for the JB potential. The 
spin-orbit splitting of the single-particle levels of A in the Op states has been 
given by 1.41, 1.57 and 1.66 for the JA, JB and NSC potentials, respectively. 
These values are very large compared with the results in Ref. [19]. It will be 
of high interest for the single-particle nature of A in l JO to be established 
experimentally. 

In the application of the UMOA some problems have still remained. The 
unitarily transformed hamiltonian [22] contains originally three-or-more-body 
interactions. The effects of these many-body interactions have been evaluated 
in the previous works of calculating nuclear properties, and we have concluded 
that the contributions of these many-body terms are much smaller than those 
of two-body interactions but sizable in some cases. We therefore should esti- 
mate the effects of these many-body interactions also in hypernuclei. 

Another refinement to be made is that, in the present calculation, the single- 
particle energies of E in intermediate states are neglected. However, the single- 
particle spectrum of E has a possibility of giving rise to some contributions 
to the AN effective interaction. Although the properties of the single-particle 
energies of E have not always been made clear, we should investigate how 
much the spectrum of E affects in the determination of the AN effective 
interaction. This problem is also a remaining task in the UMOA. 
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